A note on the resonant frequencies of rapidly rotating black holes 
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I discuss the range of validity of Detweiler's formula for the resonant frequencies of rapidly rotating 
Kerr black holes. While his formula is correct for extremal black holes, it has also been commonly 
accepted that it describes very well the resonant frequencies of near extremal black holes, and that 
therefore there is a large number of modes clustering on the real axis as the black hole becomes 
extremal. I will show that this last statement is not only incorrect, but that it also does not 
follow from Detweiler's formula, provided it is handled with due care. It turns out that only the 
first n « — log (r+ — r_)/r+ modes are well described by that formula, which translates, for any 
astrophysical black hole, into one or two modes only. All existing numerical data gives further 
support to this claim. I also discuss some implications of this result for recent investigations on the 
late-time dynamics of rapidly rotating black holes. 

PACS numbers: 



I. INTRODUCTION 

Black hole spacetimes provide an interesting arena for 
classical wave propagation: they absorb, scatter and 
sometimes even amplify incident waves. The detailed re- 
sponse of a black hole to a given incident wave depends 
on the parameters characterizing the black hole, and thus 
detection of a black hole (by, for example, gravitational 
wave detectors) requires a deep knowledge of what pre- 
cisely happens when a wave hits a black hole. 

The first numerical studies on this [l| revealed that the 
response of a black hole to a given incident wavepackct 
is dominated by an exponentially decaying ringing 
phase. What's more, the ringing frequency and decaying 
timescale of this signal depends only on the black hole 
parameters (mass, charge and angular momentum), and 
not on the particular details of the wave packet. Thus, 
a black hole, like so many other objects, has characteris- 
tic oscillation modes, called quasinormal modes (QNMs) 
(the associated characteristic frequencies are called quasi- 
normal frequencies, QN frequencies or loqn). It turns 
out that these modes are excited by almost any process 
taking place in the vicinities of a black hole, and thus a 
direct detection of black holes may well be done by iden- 
tifying this characteristic ringing |2| • This possibility has 
stimulated a great deal of research on the computation 
of QNMs and QN frequencies of black holes (see also 
Q for another recent motivation for studying QNMs in 
black hole spacetimes). The situation for Schwarzschild 
black holes is well understood: in Schwarzschild coordi- 
nates, assume that the wave has the dependence 



*(t, r, 9, 0) = e^ t+m ^Y lm (0)R(r) . 



(1) 



where Yi m {9) are just the usual spherical harmonics. I 
assume ^ is a scalar wave, otherwise the separation of the 
angular variables would have to be done with vectorial 
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or tensorial spherical harmonics. Since the Schwarzschild 
geometry is spherically symmetric, the azimuthal num- 
ber m won't play any role. The QN frequencies are those 
frequencies u> for which the field is outgoing at infinity 
and ingoing near the horizon [3. Now, it turns out that 
for each angular number I there is a discrete, infinite 
number of frequencies satisfying these boundary condi- 
tions, and that they all have a negative imaginary part, 
i.e., if we write loqn = ojr + iu>i then Wj < 0. Ac- 
cording to the time dependence in JJJ the Schwarzschild 
geometry is then stable, since any perturbation will de- 
cay exponentially with time. The QN frequencies are 
arranged by growing magnitude of imaginary part. For 
instance, considering gravitational waves, the fundamen- 
tal mode (n = 0) of quadrupole waves (I = 2) is 
ujM = 0.37367-0.08896, while the first overtone (n = 1) 
is uM = 0.34671 - 0.27391. The fundamental mode 
will in general control the signal, because it is the least 
damped, and it's quality factor Q ~ loji/luj is of order 
unity. Thus, we see that Schwarzschild black holes are 
very poor resonators, with small Q's. 

The QNM spectrum for Kerr black holes is at the 
moment well understood, for small rotation at least pj. 
Since this geometry is no longer spherically symmetric, 
the results depend on the azimuthal number m. The 
numerical studies indicate that for m > the quality 
factor increases dramatically as one increases the rota- 
tion because the magnitude of the imaginary part of the 
QN frequency decreases. Indeed, it looks like uij — > 
as a — > M, where a is the angular momentum per unit 
mass and I recall that a is bounded from above by M; 
for a = M we say the black hole is extremal. Now, what 
can we gain from this? The first obvious answer is that 
if these modes get less damped, they should be easier to 
detect. However, this is so provided they are excited to 
a measurable amplitude (sec 6] for a discussion on this 
point), which does not happen. Still, there is an analyti- 
cal result by Detweiler |jj indicating that for very rapidly 
rotating black holes there seems to be an infinity of QN 
frequencies accumulating on the real axis. Thus, it could 
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be that even though a single QNM is hard to excite, all 
of them (which cluster on the real axis, and are therefore 
nearly undamped) could conspire to give a measurable 
effect. This was investigated by Glampcdakis and An- 
dersson Q recently, who showed that this effect is indeed 
there for extremal black holes. They also suggest that 
nearly extremal black holes could also display a similar 
behavior. 

Here, I will discuss Detweiler's formula and its range of 
validity, to see whether or not this or other phenomena 
involving rapidly rotating black holes could have any im- 
portance for astrophysical black holes. Contrary to com- 
mon expectations, I find that this formula is of limited 
use for any astrophysical black hole. Indeed, it describes 
accurately only the first n < — log (r+ — r_ over- 
tones of rapidly rotating black holes. For extremal black 
holes it should describe all the modes, and therefore there 
really is an infinite number of QN frequencies clustering 
on the real axis. For non-extremal black holes, say for 
example (r + — r_)/r + ~ 10~ 10 (which is still impossible 
to find in any astrophysical scenario Q ) it describes well 
only a number n << 20 modes. 



II. THE RESONANT FREQUENCIES OF NEAR 
EXTREMAL KERR BLACK HOLES 



A Kerr black hole is characterized by two parameters, 
the mass M and the angular momentum per unit mass 
a. It has an event horizon at r+ = M + V M 2 — a? and 
a Cauchy horizon at r_ = M — V M 2 — a 2 . 

To analyze the QNMs and QN frequencies of black 
holes, one studies the wave equation that results from a 
linearization procedure. For Kerr black holes, the equa- 
tions describing linearized waves were derived by Teukol- 
sky 0- I shall not reproduce his equations here but 
merely state the results, and refer the reader to the origi- 
nal works [9j, llOj . The equations describing the evolution 
of scalar (s = 0), electromagnetic (s = 1) and gravita- 
tional (s = 2) perturbations are described by the follow- 
ing coupled equations: 



x(x + a) 



d Ri„ 
~dx r 



-{2iCjx 2 + 2x(2iuo - s - 1) - (s + l)cr + 4zr} 
-{2(2s + l)iu(x + 1) + \}Ri m = , 



dRi„ 
dx 



(2) 



1 d . . n dS . ,22 „ 2 m2 

(sin 61 — ) + ia 2 u 2 cos9 2 ~ 

sin0d0 v d6' V sin(9 2 



2ms cos 9 



sin ( 



s 2 cot 9 2 + E - s 2 )S = 0. 



2aujs cos 9 
(3) 



Here I have defined the following variables 



(7 


r + — r_ 
r+ 




{■v 




a 

2M r + 




(6) 


T 


= M (d — muj_ 


h) , 


(7) 




— ojr + 




(8) 


S 2 


= 4uj 2 - 1/4- 


A , 


(9) 



where A = E + a 2 uj 2 — 2mau> — s(s + 1). 

These equations must be supplemented by "physically 
acceptable" boundary conditions. For a scattering prob- 
lem, we allow in- and out-going waves at the asymptotic 
region (r — » oo), but only in-going waves near the hori- 
zon, i.e., the field satisfies 



R" 



1 

Z out 



-1 g2icjr* 



-oo 



(10) 



For a study of resonant frequencies, we do not want waves 
incident from infinity, and thus it must be that Z m = 0. 
I shall also be interested in the double limit a — > M, 
lu — > 77icj + , which corresponds to er — > 0, t — > 0. A 
uniformly valid solution for © in this double limit was 
first found by Teukolsky and Press following ear- 
lier work by Starobinsky and Starobinsky and Churilov 
[Tl|. Detweiler then used that solution to find the 
QN frequencies of near extremal black holes. I shall now 
re-derive the condition for the QN frequencies of near ex- 
tremal Kerr black holes, making explicit the assumptions 
that go with it. 

Let us first consider equation J5J in the limit when 
x >> max (cr, t), i.e., for large radii. Then J2J is well 
approximated by 



, d 2 Rin 



{2iu)x 2 + 2x(2iu - (s + 1))} 



dRln 

dx 



~{2(2s + l)%Q(x + 1) + \}Ri m = . 



(11) 



A solution to (|llfl that satisfies (Tp} can be written in 
terms of confluent hypergeometric functions [T^ 

R lm = Ax- s - 1/2+2lCj+lS M(l/2 + s + 2iLb + i5 1 l + 2i5,2iLbx) 
+ B(6^-6) (12) 

where A,B are constants and the notation (<5 — > —5) 
means "replace S by — S in the preceding term". 

We next turn to the case when x « 1, i.e., try to find 
a solution that is valid close to the black hole's horizon. 
In this regime, equation can be written as 



x(x + a) 



d 2 Rh 



dx 2 

-{2(2s + + X}R lt 



{2x(2iLO - (s + 1)) - (s + l)er + Air)} 
. 



R 



! in 



dx 
(13) 



This is the hypergeometric equation, and one solution 
can be written as 



r — r + 



(4) 



Rim = 2F 1 (l/2+s-2iLu+iS, \/2-2iw+s-i6, l+s-4iT/a, -x/a) 

(14) 
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It is straightforward to verify that 2-F1 — ► 1 a s x — > 0, 
which means that this solution has the desired "purely 
ingoing wave" behavior close to the event horizon. 

The key point now is that there is an overlapping re- 
gion, where both solutions are valid. Indeed, in the over- 
lap region max(cr,T) <C x <C 1 the solutions l|12|l and Ijl4(l 
both describe the solution, and can therefore be matched. 
To do the matching, take first the x — > limit of 1|12|) 
which yields, 

R lm -> Ax- s -V 2+2i&+i5 + B(S -> -<5) (15) 

Similarly, for x — > 00, l|14|) becomes, 

C((5)a;- s ^ 1 / 2+2 ^ +t ' 5 + (7(-5)x- s - 1/2+2i "- i<5 , 

(16) 

where 



C(£) = 



r(i 



4ir/cr)r(2i(5)c7 



i+l/2-2i 



r(s + 1/2 - 2iu + iS)T(l/2 + 2iuj + id + Air /a) 

(17) 

We can extract A and B by matching the solutions H15|) 
and lflU|) . 



.4 
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^ 1X1/2 + 2^ + ^ + 4^/(7) 
CW r(l/2 + ^ + z5-4zr/a) 
r(l/2 + 2»^ + »^ + 4»r/a) 
1 J r(l/2 + 2iw + i,5-4iT/a) 1 J 



On the other hand, approximating l|12l) for x — > cxd we 
get for the amplitudes Z in ,Z out , 



A- 



T(l + 2i5) 



T(l/2- s-2iuj + i5) 
B(S^-8) 



-l/2-s— 2iu-iS 



(20) 



A r(i + gg) 

r(l/2 + s + 2iw + i<5) 



(2iw) 



-l/2+s+2iti-i<5 



(21) 



To find the resonant frequencies we now impose Z in = 
0. Using iQ - CSt wc nn d that f° r these frequencies the 
following condition must be satisfied 



T(2iS)T(l + 2iS) 
T(-2i5)T(l - 2iS) 

■- (— 2ioi)cr) 



r(s + 1/2 - 2iuj - iS) 



1 2 



2l5 r(i/2- 



r( s 

2iu) - 



■l/2-2zo) + i(5) 
i<5 — Air/a) 



T(l/2 + 2iLu-iS- Air/a) 



Detweiler's formula 



(22) 



Up to this point, the only assumptions that entered the 
derivation of l|22|l were that the black hole is very rapidly 
spinning and that the lo ~ moj + . Let us now proceed 
to derive Detweiler's formula for the QN frequencies of 
rapidly rotating black holes. 



TABLE I: Values of the constant 8 2 for some (I, m) pairs, 
for au; = m/2, i.e., the extreme limit, and for scalar pertur- 
bations. Notice two important things: first, for I = m = 1 
the quantity 5 2 is negative, and thus S is purely imaginary. 
Second, for I = 6 , m = 5, 8 is purely real. Thus, although the 
rule is that for I — m the quantity 8 is real, while for other 
/ , m values it is imaginary, there are exceptions. 



/ 


m 


S 2 


I 


m 


S 2 


1 


1 


-0.4497 


3 


2 


-4.9144 


2 


2 


0.8948 


4 


4 


8.1233 


2 


1 


-4.3926 


4 


3 


-3.8792 


3 


3 


3.7552 


6 


5 


2.7752 



The quantity S is in general a complex quantity. How- 
ever, the numerical studies show that for to ~ ma+, and 
for some I, m values it is almost purely real. I would like 
to emphasize that this fact is seen numerically, it hasn't 
been analytically proved. Moreover, it seems to be true 
for most / = m modes but not only: S is real also for 
some other values of /, m, and so there is nothing special 
about I = m modes, contrary to what is stated in p. To 
be more precise, I show in Table H] the values of S 2 for 
some I, m pairs, assuming one is exactly in the extreme 
case, i.e., auj = m/2. The Table refers to the scalar case 
s = 0, although the trend is similar for other s values. 
Notice that for I = m = 1 the quantity S 2 is negative, 
and notice also that for I = 6 , m = 5 it is positive. So, 
there really is nothing special about I — m modes. 

With this cleared, I shall go on assuming that 8 is a 
real quantity. 

The left-hand side of equation 1)221) has a well-defined 
limit as a — > M and lu — > muj + . We represent that limit 
by 

LHS = qe lx . (23) 

Now, we cannot have a consistent solution unless 
t/<t — > 00 as a — > M; this is because if r/u had some 
finite limit, then so would the Gamma functions in the 
RHS of equation l|22|) . And then we would have u> ~ — , 
with X a well defined quantity. But this would violate 
our assumptions that u) ~ nuj + . Therefore, let us pro- 
ceed assuming that t/ct — > 00. Using Stirling's formula 

m, 

T[b + az] - V2^e- az (az) b+az - 1/2 , 
the right-hand side of the 1221) can be written 

RHS = (-8ujt) 2iS . 
In other words, a QNM must be a solution to 

- . 



(24) 



(25) 



f(u) = {-Slut) 2iS - qe lx = . (26) 

Using — 8ujt — pe 1 ^ we see that solutions follow from 
(remembering that 5 is real and lo — > lo + ) 



p = exp 



26 



(27) 
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C = -^iogg, 



(28) 



III. DISCUSSION AND IMPLICATIONS 



where k is any integer number. This is Detweiler's for- 
mula |2J. Since we are always assuming uj w mw+, 
this can also be written as G us t substitute — &ujt ~ 
-Am(uj - m/2)) 

wM « ^ - -L e (X-2fe 7 r)/25 cog c _ * (,c-2fcr)/M sin ( _ 

2 4m 4m 

(29) 

It is possible to prove [f| that sin£ > 0, i.e. these QNMs 
are all damped. 



B. Discussion of Detweiler's formula 

I now turn to discuss what is the range of validity of 
formula (|29|l . We have been assuming that r/er — > 00, 
but there is a number k in Ij29(l at which this condition 
breaks down. In fact, t » a implies, according to l|29|l 
that in order of magnitude one should have 



-k r + ~ r - 
e k » — . 

r+ 

Thus, there is an upper bound on k given by 

r _ r _ 



k « - log- 



(30) 



(31) 



Notice that formula is basically a lower bound for 
the fundamental mode. In fact, with we get, from 
equation {2§l that ujjM >> (r + - r_)/r+. 

This is a very restrictive bound. For instance, if we 
have (r_|_ — r_)/r+ ~ 10~ 10 , formula (|29|l can only de- 
scribe well the first n << 20 modes, where n now labels 
the overtone number. We know that accretion cannot 
spin a black hole beyond a — 0.998M ||, and so, for any 
astrophysical black hole, Detweiler's formula should only 
describe accurately one or two modes. 

Extremal black holes are accurately described by i|29|) 
for all n. In fact, all one has to do to prove this is to put 
a = in Teukolsky's equations. In this case equation 
(|14fl would be replaced by a confluent hypergeometric, 
and the final outcome for the QN frequencies would still 
be equation l|29|) . Thus, for extremal black holes, there 
is indeed an infinity of modes clustering on the real axis. 



I have re-derived Detweiler's result for the quasi- 
normal frequencies of rapidly rotating black holes, and I 
have shown to what extent it can be used. I have proved 
that for any astrophysical black hole, Detweiler's approx- 
imation is likely to be accurate for just the fundamental 
mode. In Q it was shown that the signal from extremal 
black holes goes like t~ x at late times instead of the usual 
power-law t~ 2l ~ 3 (for I = m, but I refer the reader to 
. Although they prove this for I = m modes there should 
be other modes decaying like this, according to the dis- 
cussion in this paper. Now, the main ingredient in their 
result was the fact that for extremal black holes there is 
an infinity of long-lived modes near the real axis. What 
hope is there to have this same phenomenon for non- 
extremal, but rapidly rotating black holes? As I have 
shown, this effect will not be present for any astrophys- 
ical realistic black hole, not according to the analytical 
approximation, and not according to the numerical re- 
sults in 5j. To conclude, I would like to draw attention 
to the fact that Detweiler's approximation holds good for 
any (/ , m) pair for which the quantity 5 is real. This in- 
cludes most of the I — m modes, and some other modes, 
but not all. However, from the numerical work in [j| it 
seems that all modes having positive m tend to cluster on 
the real axis, as the black hole approaches extremality. 
One is thus led to suspect that there should be some gen- 
eral argument, independent of <5, showing that all modes 
cluster on the real axis. 
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